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Abstract: We use the results from a recent investigation of hard parton-parton gravi- 
tational scattering in the ADD scenario to make semi-quantitative predictions for a few 
standard high-i?j_ jet observables at the LHC. By implementing these gravitational scatter- 
ing results in the PYTHIA event generator and combining it with the CHARYBDIS generator 
for black holes, we investigate the effects of large extra dimensions and find that, depending 
on the width of the brane, the relative importance of gravitational scattering and black 
hole production may change significantly. For the cases where gravitational scatterings 
are important we discuss how to distinguish gravitational scattering from standard QCD 
partonic scatterings. In particular we point out that the universal colorlessness of elastic 
gravitational scattering implies fewer particles between the hard jets, and that this can 
be used in order to distinguish an increased jet activity induced by gravitational scatter- 
ing from an increased jet activity induced by eg. super-symmetric extensions where the 
interaction is colorful. 
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1. Introduction 

The most exotic, and by far most discussed collider signal of large extra dimensions in the 
ADD scenario [1-3] is the copiously produced extra dimensional black holes [4-7]. While 
these are expected to come with large cross sections and characteristic signals at the LHC 
for a "natural" Planck scale of around 1 TeV, both cross section and signals suffer severely 
from uncertainties associated with quantum gravity. This provides a wonderful chance to 
probe quantum gravity, but from the point of view of verifying the scenario it is not ideal. 
It is therfore worth looking for processes which involve fewer uncertainties than decaying 
black holes. 

Other important ADD processes involve Kaluza-Klein modes, either the production of 
real ones, or the exchange of virtual ones in gravitational scattering of hard partons. The 
signal for the former involves a large missing transverse momentum and will be difficult 
to distinguish from eg. the production of stable super-symmetric particles in some SUSY 
extensions of the standard model. The later will show up as an increase of the jet cross 
sections at high energies and may, if this increase is is small, be difficult to distinguish from 
other beyond-the-standard-model effects. 

In a previous paper [8] we investigated an alternative signal for the ADD scenario, 
namely the disappearance of the high-E 1 ^ jet cross section due to the formation of black 
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holes. However, in that paper we neglected the contribution from hard gravitational scat- 
tering. 

Lately a coherent picture of gravitational scattering in the ADD model, at both low 
and high energies, was presented in [9]. In this paper we investigate the phenomenological 
consequences. Again, we will concentrate on standard jet observables to see how they 
are affected by the existence of large extra dimensions, using different choices of the model 
parameters. We will try to give a complete semi-quantitative description of the observables 
ranging from the region of perturbative gravitational scattering in the low-energy end to 
the domain of classical (non-quantum gravitational) black holes for energies above the 
Planck scale. 

While the LHC should easily discover large extra dimensions for the most natural 
choices of Planck masses and number of extra dimensions, we find situations where no 
black holes are formed and the only gravitational scattering signal could be a slight increase 
of the E± and di-jet cross section at high energies. We therefore discuss the possibility 
of distinguish such scatterings from standard QCD events by studying the different color 
topologies involved. We also suggest that such a procedure could be used at the Tevatron 
to see if an increase of the high-E± jet cross section there could be the result of the onset 
of subplanckian gravitational scattering. 

This paper is organized as follows. After a brief introduction to the ADD scenario in 
section ||, we summarize in section || the description of gravitational scatterings developed 
in [9]. We then go on to discuss the production and decay of black holes in section |I| and, 
in section [5], we present our results before discussing our conclusions in section pi 

2. Basics of ADD 

The so called ADD scenario, invented in 1998 by Arkani-Hamed, Dvali and Dimopoulos 
[1,2], aims at explaining the hierarchy problem, i.e. why the observed Planck scale at 
10 19 GeV is so large compared to the masses of the standard model particles. This is 
done by introducing a number, n, of extra dimensions in which only gravity is allowed to 
propagate. 

In order to explain why these dimensions have not yet been observed, it is assumed 
that they are compactified with some (common) compactification radius R 1 and that no 
gauge fields are allowed to propagate in the extra dimensions. Gravity, on the other hand 
is, and this renders the form of Newton's law at distances, r, much smaller than the 
compactification radius 

V(r) = S n F(n) 1 
mim 2 M™ +2 {2n) n r n+l 

Here Mp is the fundamental Planck scale, S n = 2ir n / 2 /T(n/2) is the surface of a unit sphere 
in n dimensions and T(n) is the Euler Gamma function. At distances large compared to 

We use R to denote the compactification radius rather than the compactification circumference (see 
the appendix for a discussion on conventions). 
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(2.1) 



the compactification radius we must recover the normal 3+1-dimensional form of Newton's 
law. 

v{r) - r 1 (22) 

m\m2 r 

Expressing Newton's constant in terms of the observed (3+l)-dimensional Planck scale, 
Gjv(4) ~ 1/Mp 4 , then gives the relation Mp 4 ~ Mp +2 R n between the fundamental Planck 
scale Mp and the observed 4-dimensional Planck scale, which explains how the fundamental 
Planck scale could be (almost) of the same order as the weak scale, but the observed effective 
Planck scale, Mp4, many orders of magnitude larger. 

However, this also implies that gravity should be very strong at small distances which 
opens up for the possibility of observing gravitational scattering and black holes at collider 
experiments. 



3. Gravitational scattering in ADD scenario 

Although the field theory of gravity is ultimately divergent also in more than 4 dimensions, 
an effective low-energy theory can be constructed by a perturbative treatment of the metric 
in the limit where the metric perturbation is small. A Lagrangian can be derived and 
Feynman diagrams can be constructed from it. This is done in [10, 11]. Since the extra 
dimensions are compactified, momentum occurs in each direction as multiples of some 
ground frequency, ie. as Kaluza-Klein modes. 

In a gravitational event an outgoing Kaluza-Klein (KK) mode will have some (quan- 
tized) momentum in the extra dimensions, which enters in the (3+l)-dimensional La- 
grangian as a mass term. But the KK modes can also occur as intermediate states in 
which case they have to be properly summed or, taking the continuum limit, integrated 
over. This gives rise to the integral 

mj I 

Here I enumerates the allowed momenta, mj, in the extra dimensions, m is the absolute 
value of mj, and k 2 is the momentum squared of the 3 + 1-dimensional part of the prop- 
agator. Note that this sum over KK states does imply momentum non-conservation for 
momenta transverse to the brane where the standard-model fields live, but this is not a 
complete surprise since translational invariance is broken in the bulk by the presence of 
the brane. 

3.1 Dealing with divergences 

What is worrying though, is that the field theory seems to contain divergences already at 
the tree level. However, the divergences dissappear when imposing the requirement that the 
standard model particles live on a brane, either directly by assuming a narrow distribution 
of the standard model fields into the extra dimensions [9, 12], or by a introducing a "brane 
tension" [13,14]. Both these methods gives physical effective cut-offs for the momentum 
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(mass) of the KK modes. For example, a Gaussian extension e m2 /( 2M s) Q f the standard 
model field densities into the bulk gives an "effective propagator" [9] 

D(k 2 ) = R n S n I d ™ mn ~* -rnVM* (3 2) 

J «r — m 1 

for the exchange of KK modes with four-momentum exchange k 2 . (This object, D(k 2 ), is 
here sloppily called a propagator, despite the fact that the multiplicative Lorentz structure 
is not taken into account.) 

For momentum exchange small compared to M s , the standard model momentum k in 
the propagator is irrelevant (for most m in the integral), such that s-, t-, and u-channels 
are equally efficient and the scattering is fairly isotropic. 

For y/jp* 3> M s on the other hand, the interaction is dominated by forward scattering 
via the t-channel, and an all-order eikonal calculation is necessary to ensure unitarity [9,15, 
16]. The stage is therefore set by three energy scales, the fundamental Planck mass, Mp, 
the inverse brane width (brane tension), M s , and the rest mass of the partonic scattering, 
y/s, and the phenomenology depend on their relative magnitude. It is illuminating to fix 
one of these scales an study the different kinematical regions in the plan spanned by the 
other two. This is done in fig. |l] where the (y/s, M s )-plane is plotted for Mp fixed to 1 TeV. 

Below we will successively describe the contribution from the t-, u-, and s-channels 
and the various regions in fig. |l|. 

3.2 t-channel 

As argued in the above section, we expect t-channel contributions to dominate at energies 
high compared to M s . Unitarity constraints does, however, imply that the Born approx- 
imation can not be valid for sufficiently high energies. In fact, as is argued in [15, 17], 
a completely new phenomenon occurs for scattering in more than 3 spatial dimensions; 
namely the emergence of a length scale associated with the transition from the classical to 
the quantum domain. 

Intuitively this can be understood by considering the ratios 

A ^ A M to * 

— and — (3.3) 

where 6 is the scattering angle and b the impact parameter in a scattering experiment. In 
the classical domain these ratios are both much smaller than 1. Requiring the opposite, 
and approximating 

A „ Ao h b dV(b) . 

A0 - andfl =■ ^ (3.4) 

Mv MvAb Mv 2 db y ' 

for a non-relativistic particle with speed, v, mass, M, and transverse momentum, Aq, 
moving in a potential, V, one finds, for a Coulomb-like potential V(b) = a/b, the condi- 
tion a < hv. For coupling constants close to 1, this basically implies that the relativis- 
tic and quantum mechanical regions coincide. For a more general potential of the form 
V(b) = a/b n+1 , assuming n to be positive (this is what Gausses law gives in 3+n spatial 
dimensions), the separation of the classical and quantum domain depends on the impact 
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Figure 1: The (i/i, M s ) -plane for n = 4 and Mp = 1. The straight line separating region 1 and 
4 is t/s = M s while straight line separating region 4 and 5 is the line where the real and imaginary 
parts in cq. ( |3.14| ) have equal magnitude. The power-like solid curve separating region 1 and 2 is 
yfsZ from eq. ( |3.10| ) as a function of M s and the line separating the regions 4 and 5 from region 3 
is the line where I^Bor-n^l = 1, see eq. ( 3.14 ). In the regions 1 and 2 ^/s is larger than M s , and, 
at least for ^Js 3> M s , the eikonal approximation is correct. In region 1 the eikonal is, depending 
on b, either large compared to 1 or given by eq. ( |3.8| ). In region 2 on the other hand the 6-range 
where |x| is small includes a region where it is described by cq. ( |3.9| ). In region 3 the correction 
corresponding to higher order loops is small, but in region 4 it is important and helps assuring 
unitarity. The dashed line indicates the minimal y/s (for a given M s ) at which the black hole radius 



cq. (4.1) is larger than the brane width. The plot visibly very similar for n=6. 
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Figure 2: (a) The one loop contribution corresponding to exchange of two KK modes. The 
KK modes are drawn as thick lines and standard model particles as thin lines, (b) The two-loop 
contribution. 



parameter, such that, the transition occurs at b c ~ [a/hv\^ l / n \ For gravitational coupling 
with a = Gi +n MM, this corresponds to b c ~ \G^ n M 2 / \vh)] l l n [15]. Scattering in the 
potential eq. ( |2.1[ ) is therefore expected to be mainly classical only if the impact parameter 
b is smaller than b c . 

For Mp ~ 1 TeV the parameters of this equation are such that we will see a transition 
between the classical and quantum domain at LHC, and a more careful calculation, sum- 
ming up amplitudes from ladders of t-channel exchange in fig. || to all orders is necessary. 
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This calculation was performed in [15] by simply ignoring the divergences corresponding 
to local contributions in eq. (|3,1|), and recently in [9] by a more careful analysis using the 
effective propagator eq. (|3.2[) . 

The parameter b c also corresponds to the impact parameter where the eikonal scatter- 
ing phase 



^ = h 



d 2 k ± 



ik i b i 



^4-Born( — ^J_) 



(3.5) 



becomes large compared to h. This makes perfect sense, as b c represents the impact 
parameter separating quantum mechanical and classical scattering. 

At least in the eikonal region, i.e. for small scattering angles, where there is no spin 
dependence, the Born amplitude can be written [9] 



A Bo m(k 2 = t) 



2 n-3 7r n-l M n+2 n 



dm m 



n-l 



k 2 



-m 2 /M 2 



m 



(3.6) 



where the suppression factor e m ' 2 / M s comes from implementing the requirement that the 
standard model particles live on a finite brane [9]. The same effect can be obtained by 



assuming a finite brane tension [13]. Computing the integrals in eq. (3.5) [9] then gives the 
result 



X(b) 



(20F)«M 



M 2 b 2 ^ 
— T~ ) 



(3.7) 



where the [/-functions are confluent hyper-geometric functions of the second kind. 

In the limit of large third argument, M s b 3> 1 in U, ie. impact parameters much larger 
than the brane width, % can be written 



X(b) 



br 



for b r = 



1 



sT{n/2) 
M£+ 2 



l/n 



(3.8) 



At least if b c S> 1 /M s the eikonal eq. (^^) reaches 1 in the region where it is determined by 
eq. (3.8) and b c is indeed the parameter associated with the transition from the quantum 



mechanical to the classical region. For M s small compared to Mp the brane width is more 
important and there is an energy range where the impact parameter for which \x\ reaches 
1, is given by the small argument limit in U, rather than the large argument limit, 



2s 



(20F) n M2 



Ms_ 
M P 



(3.9) 



where tp(^) is the digamma function. The transition between, \x(b)\ ~ 1 described by 
eq. ( |3.9| ), and \x(b)\ ~ 1 described by eq. (ft.8j), occurs roughly at the energy where b c = 
1/M S , and the phenomenology will therefore differ in the regions b c > \/M s and b c < 1/M S . 
Solving b c = 1/M S we find 



M n+2 n n/2 



(3.10) 
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(a) (b) 

Figure 3: When the exchanged momentum is small compared to M s , the KK propagator are 
effectively replaced by vertex factors. The diagrams in fig. ^ can then be drawn as above with only 
standard model particle lines. 

this is the line separating region 1 and 2 in fig. |l|. For t 3> M 2 the all order eikonal 
amplitude is given by 

A cik {k 2 ) = -2i8 J d 2 b ± e ik ^{e ix - 1). (3.11) 
When \x\ is large compared to 1 which, for b c S> 1/M S , happens for b < b c , the exponen- 



tiation in eq. ( 3.11 ) is important while for larger b the eikonal amplitude is approximated 
by the Born term. 

For b c < 1/M S , region 2 in fig. |l|, |x| is smaller than 1 except for very small impact 
parameters, 

b< W s eX V{ 2iM^J' (3 - 12) 

found by ignoring the digamma function in eq. (3.9). In the whole of region 1 and 2 for 



t S> Mg the gravitational cross section is obtained from the all order eikonal amplitude in 
eq. ( |3.1l| ) (although higher order corrections are only important region 1). It is given by 

If, on the other hand, -^/s <C M s , such that y/—i necessarily is small compared to M s , 
the Born amplitude is (apart from large angle spin dependences) fairly isotropic. The 
ladder-type diagrams in fig. ^ will effectively turn into </> 4 interactions as in fig. ^. 

Since the coupling grows with energy, higher order corrections will for some s become 
necessary to ensure unitarity. Summing up all contributions of the type in fig. [3|, neglecting 
large angle spin dependence, a geometric series is found [9] which helps unitarizing the cross 
section. For the 1-loop contribution we have, with P = pi + P2 as in fig. |3|, 



WMP ~ 2j q<Ms (2*r B0Tn {P/2- q f{P/2 + qf' 

1 M 2 

= A 2 Bmn ■ X with X « —(In -± + in) (3.14) 
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and higher loop corrections give similar results. Summing all ladders we obtain 

^-ladders = ^Born (1 + ^Born^ + (^Born^) 2 + ...) = ■ (3.15) 

t — ABorn^- 

Thus loop contributions of this type helps unitarizing the cross section in region 4 in fig. [I], 
defined to be the region where y/s < M s , but on-shell intermediate states in fig. [3| dominate. 
While we can not prove that these are the most important contributions, it seems likely 
as long as the cross section is dominated by on-shell states. In region 5 in fig. | this is no 
longer true since here the imaginary part in eq. ( |3.14 ) is smaller than the real part. The 



simulations performed in this paper are in the phase space region 1, 2, 3 and 4 in fig. |, 
and we use eq. ( 3.14[ ) to unitarize the cross section in regions 3 and 4 (although it's not 



important in region 3). 

3.3 u- channel 

In the regions 1 and 2, the u-channel contribution (here generally defined to be the case 
where the outgoing particle lines are crossed compared to the incoming) is small compared 
to forward t-channel scattering. In the regions 3 and 4, corresponding to y/s < M s , it 
is, however, of the same order of magnitude. In fact there is no difference between the u 
and t-channel in fig. ||. This implies that u-channel contributions run into problems with 
unitarity at roughly the same energy as t-channel contributions, and the result eq. ( 3.15| ) 
(again neglecting spin dependence) can be used also for the u-type ladders. 

The relevance of the u-channel ladders contribution is, however, significantly lowered 
by the fact that interaction among identical partons is suppressed at LHC. To get a handle 
on the importance of u-channel contribution, assume that only valence quarks contribute to 
the cross section. This is a reasonable assumption at sufficiently high momentum fractions 
and also, it will give an upper limit. The probability for the colliding partons to have 
identical flavor, spin and color is then approximately 1/10. 

3.4 s-channel 

For s <C Mg, the factor k 2 = s in eq. fl3.6|) is insignificant compared to most contributing 
KK masses and gives an amplitude similar to the t- and u-channels. There is, however, 
one complication. Due to the relative difference in sign between k 2 = s and m 2 in eq. ( |3.C| ) 
KK modes can be produced on shell. 

From the point of view of inclusive observables, these s-channel on-shell Kaluza-Klein 
states are, however, unimportant. The width of a single KK mode with mass m to decay 
into two standard model particles of energy m/2 is ~ m G^u\ giving lifetimes of order 
1000 seconds [11]. These KK modes will leave the detectors unseen. 

3.5 Phenomenology of low energy gravitational scattering 



As already mentioned section [3J, a gravitational scattering where the Kaluza-Klein mode 
is not in the outgoing state, comes with a momentum cut-off from the width of the brane, 
or from fluctuations of the brane. In the low-energy region, 3 in fig. [l], where the born 
approximation is applicable, a cut-off dependent amplitude can be used for describing the 
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interaction. From the point of view of perturbative gravitational scattering with internal 
KK modes only, this does not result in any extra parameters to describe the interaction. 
Instead it suffice to replace the Planck scale Mp by an effective Planck scale according to 



M eff 




(3.16) 



such that the Born amplitude, eq. 



(j3~6|), (neglecting spins) can be written 




(3.17) 



after integration over m, neglecting k 2 = t, u or s. In this kinematical region, gravitational 
scattering in the ADD model is a well behaved effective field theory depending on only one 
free parameter, M e g. The low-energy spin dependent footprint of the ADD scenario for 
any number of extra dimensions can then be written 



where a s is the strong coupling constant and k s , g(z), h(z) and f(z) are process dependent 
functions taking spin-dependence into account given in [18]. 

Gravitational scattering differ from standard-model and most beyond-standard-model 
processes in several ways. The experimentally most striking is probably that it increases 
with increasing energy. As the experimental situation stands today this is, however, badly 
over compensated by the decreasing parton distribution functions for high momentum fac- 
tions. The interaction is mediated by the large number of Kaluza-Klein modes, implying 
that the cross section will not have a single resonance structure, as opposed to cross sec- 
tion signatures of most other beyond-standard-model particles. Due to the different spin 
dependence of gravitational scattering the angular distribution will also differ. 

We will here consider another difference, namely that contrary to the main contribution 
to inclusive cross sections, both in the standard model and in super-symmetric extensions, 
the gravitational interaction is colorless. As we shall see, this implies noticeable differences 
in particle multiplicity outside the jets. 

4. Black holes in ADD scenario 

Black holes with mass large compared to the fundamental Planck scale, but with radius 
small compared to the compactification radius are expected to behave much like extra di- 
mensional versions of astronomical 3+1-dimensional black holes. The Schwarzschild radius 
is given by 




(3.18) 



?"Sch = 



i m bh sr(B±3) 



n+1 




y/irM-p M P n + 2 
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and the temperature is given by [19] 

m n+ 1 



(4.2) 



4vrr Sc h 

Note that small black holes are hotter. 

A major difference between black holes in the ADD scenario and ordinary 3- 
dimensional black holes is that ADD black holes do not radiate gauge fields into most 
of phase space, since only gravity is allowed to propagate in the extra dimensions. One 
may believe that this would lead to almost no radiation on the brane (where gauge fields 
and, hence, also we live) as the bulk phase space is much larger. However, it has been 
shown that this is not necessarily the case [4]. 

Since we are considering the non-idealized situation of a finite brane width we must 
also consider the implications of this on black hole production. In particular, a natural 
requirement is that the brane is not more extended than the black hole, leading to the 
condition rs c h < for the formation of black holes. As we will see this prevents black 

holes from appearing at the LHC for sufficiently small M s . 

On the other hand, if M s is large, we may, with increasing y/s, go directly from the Born 
region 3 in fig. [I] to black hole production. This should be worrying since the black holes 
are treated semi-classically but the gravitational scattering in region 3 is purely quantum 
mechanical. It is reasonable that the black holes should start behaving classically when 
the Compton wave length is of the same order as the black hole radius, but it would have 
been more comforting to only study black hole production in region 1 in fig. |l], where the 
gravitational scattering is mainly classical already at lower energies. This represents a 
genuine quantum gravity uncertainty. 

Already at a classical level the cross section for black hole creation is subject to sig- 
nificant uncertainties. This is basically due the fact that it does not suffice to consider 
the colliding objects, but in addition the curvature of space-time far outside the black hole 
need to be calculated. Classical numerical simulations for black hole formation in extra 
dimensions have been performed in [20,21] with the result that the geometric cross section, 
irr 2 , should be multiplied with a factor ~ 0.7 — 3, increasing with the number of extra 
dimensions. For this paper we have, however, chosen to keep the constant at 1. 

As the black holes considered here are formed from partons inside the protons there 
is also an uncertainty from the usage of parton distribution functions for an essentially 
non-perturbative process [22] . (A discussion about the effects of quantum fluctuations 
based on wave packages can be found in [23,24].) 

Then there is the question of the onset of black hole production. It can be argued 
that no black holes should be formed below (roughly) the Planck scale as the uncertainty 
principle would forbid sufficient localization of the partons. But precisely when does black 
holes begin to form? 

We consider first the condition that the black holes have to be well localized in our 
ordinary dimension. Looking at the momenta of the incoming partons in their combined 
rest frame it is reasonable to require that their wavelength, oc 2/ \fs, is less than rs c h- 
The corresponding requirement in the transverse direction gives the requirement: Aj_ oc 



- 10 - 



1/pt < r^ch- Clearly one can argue about the proportionality constant. We have chosen 



M min = 2/r Sc h(M min ). (4.3) 



Combining this with the expression for the Schwarzschild radius eq. (4.1) we get 



M min = M P 



8r(^±3) 



i 

n + 2 



(4.4) 



Numerically the value of M m ; n is then approximately twice the Planck mass. 

As we consider a finite brane width we must add the condition r^ c h > 1/M S , leading 
to the minimal mass 



Mg +2 (2 + n)^ 

8 r Ml 



M min2 = F rQ n — -tt — . (4.5) 



Again one can argue about the proportionality constant. While we take into account 
the effects of a finite brane, we do not consider the dynamics governing the brane, and 
possibly describing its width, although this may have significant effects on the spectra 
observed [25,26]. 

Once a black hole has formed it is believed to lose most of its geometric asymmetries in 
a short period referred to as the balding phase. This phase leaves a black hole whose only 
geometric asymmetry can be described by one angular momentum parameter. However, 
it turns out that this angular momentum tends to be lost rather quickly via Hawking 
radiation, such that the black hole (apart from gauge charges) can be described by the 
Schwarzschild metric. 

Neglecting the gauge charges, which in the case of electromagnetism has been shown to 
have a modest influence [27], the disappearance of the black hole would be well described 
by Hawking radiation if the black hole was much heavier than the Planck mass, and if no 
brane effects, such as the black hole recoiling of the brane [28,29], or interacting with the 
brane [25,26] is taken into account. The problem is that most collider-produced black holes 
will not be much heavier than the Planck mass. 

For a hole which is not heavy compared to the Planck mass one cannot treat the 
metric as a static background for the emitted quanta, the back-reaction of the quanta 
to the metric should be taken into account and this is not done in the derivation of the 
Hawking radiation [30]. Also, at some point, the lifetime of the black hole becomes shorter 
than its radius. This makes it difficult to talk about a thermalized black hole. 

Considering all of this, it should not come as a surprise if black holes where observed 
with spectra which differs significantly from that expected from eq. ([4.2|). 



5. Results 

We have used the amplitudes for gravitational scattering for the different regions in fig. p] 
presented above to reweight the standard QCD 2 — ► 2 scatterings in the PYTHIA (version 
6.2 [31]) event generator. In the regions 1 and 2 we have used the (elastic spin-independent 
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t-type) all order eikonal cross section from eq. (|3.5[) and eq. ( |3 .11 ). In the regions 3 and 



4 we have used the spin dependent Born amplitude [32] corresponding to eq. ( 3.18 ), and 



higher order corrections according to eq. ( |3.15 ). In region 3 the higher order corrections 



are small, but in region 4 they are essential. In the case of particle-antiparticle scattering, 
such that the scattering can be mediated via the s-channel, we have "unitarized" also the 
s-channel contribution in region 4 (and 3) using eq. (3.15), despite that fact the the s-type 



ladders, diagrams in fig. || rotated by n/2 do not have on-shell intermediate standard model 
particles. There are thus several fundamental uncertainties associated with gravitational 
scattering in region 4. First, we use the spin dependent Born amplitude, but we do not 
take spin dependence consistently into account in eq. ( |3.15 ) since we use the same A-Q oin (t) 



everywhere in all ladders. Second, we suppress s-type contributions in the same way as 
u- and t-type. (Note that we call the ladders in fig. [| and fig. || t-type, sometimes these 
diagrams are referred to as s-channel, since the resummation is in s.) 

The different treatments in the various regions means that we could expect a discontin- 
uous transition when yfs is increased, such that we cross the line yfs = M s in fig. |]. As long 
as the Born approximation is applicable (regions 2 and 3), this transition just corresponds 
to starting neglecting spin-dependence in region 2. If the transition is between region 4 
and 1, the situation is, however, worse due to fundamental uncertainties associated with 
region 4. 

For each generated 2 — > 2 scattering we also change the color flow between the scattered 
partons with a probability <tadd/(<7ADD + cqcd) to reflect the colorless nature of the 
graviton exchange. The resulting partonic state is then allowed to evolve a QCD cascade 
and is finally hadronized to produce fully simulated hadron- level events. Where relevant, 
we have also added multiple soft and semi-hard QCD scatterings to simulate the underlying 
event according to the model implemented in PYTHIA [33]. 

In addition, we have used the CHARYBDIS [34] program to simulate the production 
and decay of black holes as described in section |] and in [8]. To ensure that the energy 
is sufficiently localized, in our ordinary dimensions and in the extra dimensions, we have 
required a minimal black hole mass according to eq. (4.4) and eq. (|4.5Q . We also use the 



Schwarzschild radius to cut off any QCD and gravitational 2^2 scatterings in region 3 
and 4 for large enough masses and transverse momenta as discussed in [8]. In region 1 
(and 2) we use the impact parameter description defined via eq. ( |3.7| ) to turn off gravita- 
tional interactions at distances smaller than rg c h- Clearly this simpleminded approach of 
turning of gravitational scattering should not be seen as the final word. In particular our 
understanding of gravitational scattering, and hence its turnoff, is limited in region 4. 

We limit our investigation to two standard inclusive high-E'x jet observables [8,22,35], 
namely the -Ej_-spectrum of the highest- E± jet in an event, and the distribution in invariant 
mass, Mjj, of the two highest--Ej_ jets in an event. As high-Sj_ jets will be a part of almost 
any signal of new physics at the LHC, such observables will be measured early on after 
the start of the experiments and it is also where one would expect gravitational scatterings 
to contribute. We use a simple cone algorithm 2 with a code radius of 0.7, assuming a 



2 The GETJET algorithm originally written by Frank Paige. 
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Figure 4: (a) The -Ej_-spectrum of the highest E± jet in an event, and (b) the invariant mass 
spectrum of the two highest E± jets in an event at the LHC. In both cases M e g = 1 TeV with 
4 extra dimensions and M S /M P = 2 (M P w 0.9 TeV, M s w 1.8 TeV). The long-dashed lines 
are the contribution from QCD scatterings, short-dashed lines the contribution from gravitational 
scatterings, dotted lines the contribution from the decay of black holes and full lines the sum of all 
contributions. 



calorimeter covering the pseudo-rapidity interval, \rj\ < 2.5, and requiring a minimum E± 
of 100 GeV for the resulting jets. We have checked that our results do not depend much 
on the algorithm chosen. 

In fig. |H we show generated the E± max and Mjj 
distributions at the LHC for the case of four extra 
dimensions, Af c g = 1 TeV and M s /Mp = 2 (see 
table H for the resulting values of Mp and M s ). 
In the -Eimax spectra we see that the cross section 
is dominated by QCD scatterings at low E_\_as ex- 
pected, followed by an intermediate region where 
gravitational scattering becomes important before 
black-hole production starts dominating the cross 
section at large E±. For the Mjj-spectrum, the 
situation is different, and the gravitational scat- 
terings dominates at large masses. 

Modulo effects of the parton densities we ex- 
pect both gravitational scattering and black-hole 
production to increase with energy. For black holes 
one may naively not necessarily expect to find high-E± jets, as energetic quanta are Boltz- 
mann suppressed in the Hawking radiation. However, it turns out that the large cross 
section for a black hole to form at high s, multiplied with the small probability for the 
Black hole to radiate extremely energetic quanta, may dominate over the non-black hole 
cross section for rather large transverse momenta [8]. (Even if well localized QCD and 
gravitational scattering events are not suppressed due to black hole production.) These 
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Table 1: The different values of M e g, 
number of extra dimensions, n, and the 
ratio of M s /Mp used in the simulations 
together with the resulting approximative 
values of Mp and M s . The masses are all 
given in units of TeV. 
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Figure 5: The contribution of different regions in figure [I] to the di-jct mass spectrum from 
gravitational scatterings at the LHC with M c ff = f TeV, 4 extra dimensions and M s /Mp = 1 (a), 
2 (b) and 4 (c). In all cases the full line is the sum of all contributions and the contributions 
from regions 1, 2, 3 and 4 is given by the long-dashed, short-dashed, dotted and dash-dotted lines 
respectively. 



extremely energetic quanta do, however, not obey the semiclassical approximation in the 
Hawking radiation derivation, and are therefore associated with large uncertainties. 

For large E±, however, the gravitational scattering events, just as the QCD ones, may 
be localized inside the Schwarzschild radius and will collapse into a black hole. 

In fig. H we show the Mjj-distribution of gravitational scatterings only, divided into the 
contributions from the different regions in fig. |]. Keeping M e g = 1 TeV and the number of 
extra dimensions (4) fixed, we vary M s /Mp and find that the contribution from region 1 
dominates except in the low-mass regions below M s . The transitions between the regions 
are not sharp, mainly due to the smearing introduced by shower, hadronization and the 
jet reconstruction. This smearing hides the fact that the transition between \fs > M s 
and y/s < M s is discontinuous in the distribution of the generated s. In the case this 
transition occurs between region 2 and 3, where the Born approximation is applicable, the 
discontinuity is not even visible in the generated s-distribution. If the transition occurs 
between region 4 and 1, a discontinuity can, however, be seen. 

We note in table [t] that, although M e g is kept fixed, giving the same amount of gravita- 
tional scattering at \fs <C M s , increasing the ratio M s /Mp will increase both M s and Mj>. 
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Figure 6: The same as fig. [§a, but with M s /Mp = 1 (a), 4 (b) and 0.5 (c). 



And since black-hole production depends of Mp and M s differently via eqs. ( |4.1| ), ( fhq ) 
and (|4.5| ), we can vary the relative importance of gravitational scattering and black- hole 
production by varying M s /Mp. Hence we see in fig. |6|a that lowering M s /Mp to 1, the 
gravitational scattering will never give a sizeable contribution to the -Ej_ max -distribution, 
while increasing the ratio to 4 (fig. |6|b) results in the gravitational scattering dominating 
the cross section further out in E± as compared to fig. ^|a. In fig. ||c we decrease the ratio 
even further to 0.5 which results in a brane thickness so large that black holes can never 
be formed at the LHC, and the only indication of the presence of extra dimensions in the 
E± spectra is a slight increase in the cross section for large £i max . 

A similar effect can be obtained by increasing the effective mass, while keeping the ratio 
M s /M p fixed, hence increasing both Mp and M s . This is done in fig. 0a and, again, the 
only visible effect of the extra dimensions is from gravitational scattering in the high-£Jj_ 
region. On the other hand we see in fig. 0b how increasing the number of extra dimensions 
to 6, keeping M e g- = 1 TeV, gives a negligible contribution from gravitational scattering 
to the E'xma.x-distribution, which instead is completely dominated by the decay of black 
holes. 

If large extra dimensions exist, one would hope that the scales are such that they 
would be easily discovered at the LHC by, eg. the striking signature of a decaying black 
hole. However, it is easy to see how nature could conspire, such that the only signal in 
the E±_ spectrum would be a slight increase of the high-E'x jet cross section. There are, of 
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Figure 7: The same as fig. [I]a, but with (a) 4 extra dimensions, M e s = 4 TeV, M s /Mp = 4 and 
(b) 6 extra dimensions, M c g = 1 TeV, M s /Mp = 2. 
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Figure 8: The average number of charged particles outside the jet-cones in the central rapidity 
unit between the two hardest jets in events corresponding to fig. ||c. The full line is for all events 
while the long-dashed and short-dashed are for QCD and gravitational scatterings respectively. In 
(b) only events with a minimum pseudo-rapidity difference of one unit between the two highest E± 
jets are included, while in (a) there is no such requirement. 



course, other signals, such as the production of real gravitons, showing up as large missing 
transverse momenta. But such signals could also be the result of other possible beyond- 
the-standard-model scenarios. In any case, it would be desirable to be able to distinguish 
gravitational scatterings from standard QCD events. One obvious difference is that the 
exchange of a graviton is colorless in contrast to a QCD scattering. This will necessarily 
give rise to a different color topology in gravitational events as compared to QCD ones. In 
particular one would expect the appearance of so-called rapidity gaps between the jets in 
gravitational scattering events. Although these gaps may be filled by secondary soft and 
semi-hard QCD scatterings, one may still expect a lower activity between the jets in such 
events. 

In fig. [8|a we show the average number of charged particles with a transverse momentum 
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Figure 9: (a) The -E^-spectrum of the highest E± jet in an event, and (b) The average number of 
charged particles outside the jet-cones in the central rapidity unit between the two hardest jets at 
the Tevatron for M e g = 700 GeV, M s /Mp = 4 and 4 extra dimensions. The full line is for all event 
while the long-dashed and short-dashed are for QCD and gravitational scatterings respectively. 



above 0.5 GeV outside the jet cones in the middle unit of pseudo-rapidity between the two 
hardest jets as a function of E± max . Hence, we count only charged particles, c, with 

p ±c > 0.5 GeV, 
AR cl ,AR c2 > 0.7 and 
rji + t?2 



< 0.5, (5.1) 



where r\ c and 77, are the pseudo rapidities of the particle and (the center of) jet i respectively 
and AR c i is the distance between the particle and jet i in the pseudo-rapidity-azimuth- 
angle (rj, (p) plane. In this simulation we have included multiple interactions in PYTHIA to 
simulate the underlying event. 3 We see that the expectation from QCD events is around 
10 particles, while for gravitational scatterings the average is around 5. With sufficient 
statistics it could therefore be possible to observe the decrease in the number of charge 
particles with increasing i?j_ max as gravitational scatterings starts to dominate. In fig. ||a 
we have not required a large rapidity separation between the jets. Doing so would increase 
the effect, as shown in fig. |8]b, but on the other hand the statistics would decrease. 

We note that the absolute numbers in fig. |8| is very sensitive to the modeling of the 
underlying event, which is very difficult to predict for the LHC. The underlying event 
should, however, give the same contribution to both scattering types, and the difference 
between the two should be fairly well predicted by PYTHIA. 

At the Tevatron there was an indication of an excess of the cross section for very 
high E± jets as compared to the QCD prediction [37,38]. Although a re-evaluation of the 
uncertainties due to the parton density parameterizations used in the QCD calculations 
has brought this excess within the limits of the statistical and systematical errors, it is still 



'Using parameter settings according to the so-called Tune- A by Rick Field [36]. 
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intriguing that such an excess could be the signal of the onset of gravitational scattering 
due to the presence of large extra dimensions. 

In fig. ||a we show our prediction for the E± mSuX distribution for n = 4, M e g = 700 GeV 
and M s /M-p = 4 at the Tevatron. The parameters were chosen so that the excess above 
standard QCD production is approximately within the statistical and systematical uncer- 
tainties of the corresponding Tevatron measurement. In fig. we then show the average 
number of charged particles outside the cones (same as in fig. ||, but counting charged 
particles with transverse momenta down to 0.25 GeV). The decrease in the region where 
gravitational scatterings become important is significant, although it may be difficult to 
get enough statistics to measure it even for Run-II at the Tevatron. However, it is not com- 
pletely inconceivable that by finding a more sensitive observable of the color structure, we 
could be able to see the first indication of large extra dimensions already at the Tevatron, 
before LHC is switched on. 

6. Conclusions 

We have studied gravitational scattering and black hole production at the LHC and the 
Tevatron in the ADD scenario assuming that brane on which the standard model fields live 
have a finite width. We found that the relative importance of gravitational scattering and 
black hole production is sensitive to this width, and that for large widths the extension 
of the standard model particle fields into the bulk, may prevent black holes from forming 
since the energy may not be sufficiently localized within the black hole radius. 

A wide brane corresponds to a low cut-off for virtual Kaluza-Klein modes (the brane 
width and the KK cut-off, M s , are inversely related via a Fourier transform [9]) and there- 
fore results in weaker gravitational interaction in the non-classical regions. It is thus 
possible for nature to conspire, by choosing a low M s , such that neither much gravitational 
scattering or black holes is observed at the LHC. In this case processes involving the pro- 
duction of on-shell KK modes resulting in missing E± may become important observables. 

In our simulations we have used values of Mp, M s and the number of extra dimensions, 
n, which we believe have not yet been excluded by experiments (see eg. [39,40] for recent 
reviews). Most of these limits are only relevant to Mp but restrictions on M s could be 
obtained by considering processes involving both virtual and real KK modes. 

In the case of low M s , only a weak increase in the jet spectra could be observed at 
LHC, and the signal of missing E± could be the result of SUSY. We point out that the 
colorless nature of gravitational scattering could be a way of distinguishing gravity induced 
events from other beyond-standard-model extensions. In fact this method could be used 
to indicate if an excess of jet activity at high transverse energies at the Tevatron is a result 
of gravitational scattering. 

A. Appendix 

There are at least four definitions of the Planck mass. Often one have to understand which 
definition an author uses by the relation of the Planck mass to the 4-dimensional Newton's 
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constant Gjv(4) or to the Schwarzschild radius of a black hole. The process of hunting down 
constants is further complicated by the use of different definitions of the compactification 
radius, many authors [10] mean by the compactification radius rather the compactification 
circumference, here denoted L, whereas others really mean the radius, R = L/(2ix). In 
order to simplify comparison between the different conventions we here state the relations 
between the Planck masses Mp (used here) and in [34], used in [11,15], Mq, Ms used 
in [10] and the 4-dimensional Newton's constant Gjy(4), the relation between the Planck 
masses and Schwarzschild radius, and the relations of the Planck masses to each other. 
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